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ABSTRACT 

We  extend  and  modify  recent  efficient  techniques  for  half-plane  range 
searching,  to  obtain  algorithms  for  answering  efficiently  several  types  of 
queries  on  planar  arrangements  of  segments  and  certain  related  structures. 
For  example,  we  show  that,  given  n  segments  in  the  plane,  we  can  preprocess 
them  in  0(n  log  r)  randomized  expected  time  and  0(n)  space,  so  that,  given 
a  query  ray,  emanating  from  an  arbitrary  point  in  an  arbitrary  direction,  we 
can  determine  the  first  segment  (if  any)  hit  by  the  ray,  in  time  OCn^'^"*^'),  for 
any  €>0.  Other  results  involve  "implicit"  point  location  in  an  arrangement  of 
(possibly  intersecting)  triangles,  polygon  placement  queries,  and  certain  res- 
tricted types  of  ray  shooting  and  hidden  surface  removal  for  3-dimensional 
polyhedral  scenes. 


1.   Introdoction 

We  begin  with  the  following  somewhat  unusual  opening  remarks.  The  results  of  this 
paper  have  been  conceived  and  developed  in  the  summer  of  1987,  but  for  various  reasons 
were  left  unpublished  for  almost  two  years.  Meanwhile,  there  has  been  significant  progress 
on  the  problems  studied  here.  In  particular,  a  recent  work  by  Agarwal  [Ag]  present  tech- 
niques that  achieve  query  time  roughly  OC^{n)) ,  and  involve  deterministic,  rather  than  ran- 
domized, preprocessing  (however,  the  preprocessing  time  no  longer  remains  close  to  linear). 
Agarwal's  results  are  based  on  a  decomposition  technique  that  is  different  from  the  one  used 
in  this  paper.  Still,  many  of  Agarwal's  arguments  are  taken  from  this  paper.  As  a  public  ser- 
vice, to  make  our  results  somewhat  more  accessible,  we  have  decided  to  elevate  our  paper 
from  the  status  of  "unpublished  manuscript"  to  that  of  a  technical  report. 

Let  G  be  a  collection  of  n  segments  «i,  ...,<„  in  the  plane,  and  let  A  =  A  (G)  be  the 
arrangement  of  these  segments,  that  is  the  planar  map  whose  vertices  are  the  endpoints  and 
the  intersection  points  of  the  segments,  whose  edges  are  maximal  connected  subsegments  of 
the  segments  not  containing  any  vertex,  and  whose  faces  arc  the  connected  components  of  the 
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complement  of  the  union  of  the  segments.  In  this  paper  we  consider  a  variety  of  problems  in 
which  we  wish  to  process  such  an  arrangement  A  so  as  to  be  able  to  answer  certain  queries 
involving  A. 

Typically,  the  queries  we  are  interested  in  are  trivially  answered  in  O(log  n)  time,  if  we 
are  willing  to  calculate  and  store  A  explicitly.  Since  the  combinatorial  complexity  of  A  can  be 
quadratic  in  n  in  the  worst  case,  this  might  be  too  expensive,  and  we  therefore  seek  alterna- 
tive methods  of  more  efficient  but  implicit  representation  of  A,  which  still  facilitates  reason- 
ably efficient  query  handling,  although  not  as  efficient  as  the  (optimal)  O(logn)  time  that 
can  be  achieved  with  explicit  representation. 

As  a  typical  example,  consider  the  following  implicit  point  location  problem.  Given 
such  an  arrangement  A,  preprocess  it  so  that,  given  any  query  point  x,  we  can  calculate 
quickly  certain  properties  of  the  face  of  A  containing  x.  For  example,  suppose  A  is  formed  as 
the  overlay  of  n  (possibly  intersecting)  triangles,  and  the  property  in  question  is  how  many 
triangles  contain  the  query  point  x,  or  more  generally,  if  each  triangle  is  associated  a  certain 
weight,  what  is  the  sum  of  the  weights  of  the  triangles  covering  x,  etc.  Clearly,  by  calculating 
the  arrangement  A  induced  by  these  triangles,  and  by  further  preprocessing  it  for  efficient 
point  location,  it  is  easy  to  answer  such  a  query  in  (optimal)  0(log  n)  time.  Can  we  obtain 
faster  preprocessing  and  more  compact  representation  of  A  that  will  still  enable  us  to  answer 
these  queries  reasonably  efficiently? 

In  this  paper  we  show  that  this  is  indeed  possible.  We  extend  and  modify  several  recent 
techniques,  based  on  partition  trees  for  sets  of  points  in  the  plane,  which  have  been  originally 
developed  for  half-plane  and  polygon  range  searching  [HW],  [EW],  [DE].  Using  such  parti- 
tion trees,  we  are  able  to  solve  the  implicit  point  location  problem,  as  well  as  many  other 
related  ones,  using  0{n  polylog(n))  preprocessing  time  and  space,  and  ©(n^'^"^')  query  time, 
for  any  €>0. 

Perhaps  the  most  interesting  application  that  we  have  is  that  of  ray  shooting  in  non- 
simple  polygons,  or,  more  generally,  in  an  arbitrary  arrangement  of  (intersecting)  segments. 
In  this  problem  we  wish  to  preprocess  the  given  segments  so  that,  given  any  query  ray  p, 
emanating  from  an  arbitrary  point  in  an  arbitriiry  direction,  we  can  determine  the  first  seg- 
ment (if  any)  hit  by  p.  Using  a  somewhat  complex  partition-tree  technique,  and  an  interesting 
development  of  a  linear  ordering  of  the  segments  which  is  consistent  with  the  order  in  which 
they  arc  hit  by  certain  types  of  rays,  we  obtain  an  algorithm  that  uses  0{n  polylog(n)) 
preprocessing  time  and  space,  and  achieves  0(«^'^"^*)  query  time,  for  any  €>0.  In  the  case 
of  a  simple  polygon,  Chazelle  and  Guibas  [CG]  have  developed  a  ray  shooting  procedure  that 
uses  0(n  log  log  n)  preprocessing  time  and  0{n)  space,  and  achieves  optimal  O(log  n)  query 
time;  however,  no  efficient  solution  was  known  for  the  non-simple  case  (see  some  further 
discussion  on  this  problem  in  Section  3). 

We  also  present  additional  applications  of  our  technique,  including 

(i)  Polygon  placement  queries.  Here,  given  an  arbitrary  polygonal  object  P  and  a  collection  of 
polygonal  obstacles,  we  wish  to  determine  whether  a  specific  query  translation  of  P  will  move 
it  to  a  free  placement,  in  which  it  does  not  intersect  any  obstacle. 

(ii)  Implicit  hidden  surface  removal  in  three  dimensions .  Here,  given  certain  collections  of  tri- 
angles in  three  dimensions,  we  wish  to  preprocess  them  so  as  to  be  able  to  answer  efficiently 
ray  shooting  queries,  each  asking  for  the  first  triangle  (if  any)  hit  by  a  given  query  ray.  This 
is  a  central  subproblem  of  the  ray-tracing  problem  in  computer  graphics  [SML].    ' 


An  important  variant  of  these  problems  is  when  all  the  queries  are  known  in  advance, 
as  is  the  case  in  many  applications.  We  call  this  variant  the  "batched"  version  of  these  prob- 
lems. We  show  that  batched  problems  can  be  solved  considerably  more  efficiently  than 
unbatched  ones,  by  constructing  "customized"  partition  trees,  which  are  tailored  to  handle 
efficiently  only  the  given  queries,  disregarding  efficiency  considerations  for  any  other  query. 

An  important  theme  of  this  paper  is  to  point  out  the  versatility  of  the  partition  tree  tech- 
niques, and  their  applicability  to  problems  which  do  not  always  seem  directly  related  to  range 
searching.  In  particular  we  note  that  the  batching  techniques  introduced  in  [EGSh]  can  be 
applied  to  a  large  variety  of  problems,  much  beyond  the  specific  applications  studied  in 
[EGSh].  This  theme  has  also  been  observed  and  exploited  by  Dobkin  and  Edelsbrunner 
[DE],  where  a  weaker  form  of  partition  trees  (that  is,  "conjugation  trees")  is  used  for  space 
intersection  queries,  in  a  manner  similar  to  the  one  used  here. 

The  paper  is  organized  as  follows.  In  section  2  we  present  efficient  techniques  to  handle 
(certain  variants  of)  the  implicit  point  location  problem.  Section  3  analyzes  the  planar  ray 
shooting  problem.  Section  4  presents  techniques  for  handling  certain  three-dimensional  ray 
shooting  problems.  Section  5  handles  polygon  placement  queries.  Concluding  comments  and 
discussion  is  given  in  Section  6. 

2.    Implicit  Point  Location 

The  planar  point  location  problem  is  one  of  the  most  fundamental  problems  in  computa- 
tional geometry.  It  was  studied  in  many  papers,  culminating  in  several  optimal  algorithms 
[Ki],  [EOS],  [ST].  In  this  problem  one  is  given  a  planar  subdivision  (or  map)  M,  consisting 
of  n  faces  (and  0(n)  edges  and  vertices),  and  the  goal  is  to  preprocess  M  into  a  data  structure 
which  supports  fast  point  location  queries,  in  which,  given  a  query  point  x,  we  wish  to  deter- 
mine the  face  of  M  containing  x.  The  algorithms  in  [Ki],  [EGS],  [ST]  achieve  this  goal  in 
0{n  log  n)  preprocesing  time,  0(n)  storage,  and  O(log  n)  query  time;  if  the  given  map  has 
some  favorable  properties  (e.g.  is  a  convex  subdivision)  then  preprocessing  can  be  done  in 
linear  time  [Ki],  [EGS]. 

A  related  (and  simpler)  problem  is  when  the  points  to  be  located  in  M  are  all  known  in 
advance.  In  this  case,  if  m  points  are  to  be  located,  one  can  use  a  simpler  sweeping  algorithm 
whose  complexity  is  0((m -t-/i)  log  n)  [Pr].  This  complexity  is  comparable  with  what  one 
could  get  using  the  preprocessing  approach,  but  the  algorithm  is  much  simplified.  We  call  this 
version  of  the  problem  the  "batched"  version. 

In  this  section  we  study  the  following  generalization  of  the  problem.  Suppose  the  map 
M  is  not  given  to  us  explicitly,  but  rather  it  is  defined  as  the  arrangement  of  n  (usually  inter- 
secting) geometric  objects  of  some  simple  shape.  We  will  assume  that  the  given  objects  are 
line  segments  (or  sometimes  full  or  half  lines);  this  will  also  enable  us  to  handle  general 
polygonal  objects.  As  will  be  demonstrated  below,  this  situation  arises  in  many  applications, 
including  several  visibility  and  placement  problems.  A  naive  approach  to  this  modified  prob- 
lem is  of  course  to  calculate  the  arrangement  M  formed  by  the  given  segments,  and  then 
preprocess  it  for  point  location  as  above.  Calculating  the  arrangement  can  be  accomplished  in 
time  O(n^)  and  space  O(n^),  e.g.  by  calculating  the  arrangement  of  the  lines  containing  the 
segments  as  in  [EOS].  However,  since  the  input  size  is  only  0{n),  we  face  the  challenge  to 
improve  upon  this  naive  approach,  and  obtain  subquadratic  solutions. 

This  goal  is  not  always  possible.  For  example,  if  the  number  m  of  points  eventually  to 
be  located  in  Af  is  a  n^,  the  above  method  is  optimal,  and  no  subquadratic  performance  is 
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possible.  Even  in  this  case,  however,  there  still  remains  the  problem  of  reducing  the  space 
requirement  of  the  algorithm  to  subquadratic. 

Another  issue  at  hsmd  is  why  point  location  in  Af  is  desired.  Often  this  is  needed  for 
determining  some  "local"  property  of  the  query  point.  For  example,  if  M  is  defined  as  the 
arrangement  of  n  triangles,  then  we  might  want  to  determine,  for  a  query  point  x,  whether  it 
is  contained  in  the  union  of  these  triangles,  or  how  many  triangles  contain  x,  or,  if  each  trian- 
gle is  assigned  some  weight,  what  is  the  sum  of  the  weights  of  the  triangles  containing  x,  or 
the  minimum -weight  triangle  containing  x,  etc.  Similarly,  for  an  arrangement  of  segments, 
we  might  seek,  for  a  query  point  x,  which  segment  lies  directly  below  x.  We  think  of  these 
problems  as  being  local,  because  they  do  not  seek  precise  knowledge  about  the  position  of 
the  query  point  x  in  the  entire  arrangement,  but  rather  confine  themselves  to  combining 
pieces  of  data,  each  obtained  by  "confronting"  x  with  one  of  the  objects,  in  some  associative 
manner  (which  can  be  done  naively  in  linear  time). 

This  section  assumes  that  point  location  in  M  is  indeed  required  for  such  a  "local"  task. 
We  show  that  in  many  cases  these  problems  can  be  solved  efficiently  in  subquadratic  time 
and  space,  by  avoiding  explicit  calculation  of  M,  and  by  storing  instead  some  implicit 
representation  of  it. 

2.1.   Statement  of  the  Problem 

For  the  sake  of  exposition,  we  will  study  the  following  specific  problem.  Given  a  col- 
lection G  =  {*i,  .  .  .  ,e„}  of  n  segments  in  the  plane.  With  each  segment  «€G  we  associate  a 
function  <!>,  defined  on  the  entire  plane  which  assumes  values  in  some  associative  and  com- 
mutative semigroup  S  (whose  operation  will  be  conveniently  denoted  by  "+").  We  assume 
that  evaluation  of  4),(x)  at  a  point  x,  as  well  as  addition  of  two  elements  in  5,  takes  constant 
time.  Two  versions  of  the  problem  concern  us.  In  the  preprocessing  version,  we  want  to 
preprocess  the  data  into  a  data  structure  so  that,  given  any  query  point  x,  we  can  calculate 
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efficiently.  In  the  batched  version,  we  are  given  m  points  p,,  .  .  .  ,p„,  and  we  want  to  calcu- 
late <J)(pi) <j>(Pm)  efficiently. 

The  functions  4),  and  4>  are  assumed  to  satisfy  certain  properties,  which  make  their  cal- 
culation closely  related  to  point  location.  One  assumption  is  that  <(>  is  constant  on  each  face  of 
the  arrangement  A  formed  by  the  segments  in  G  (or  at  least  that  <})  is  constant  on  each  sub- 
face  obtained  by  some  standard  triangulation  of  other  decompositions  of  A).  We  also  assume 
that  these  constant  VEilues  of  4>  can  be  obtained  at  constant  cost  per  face  as  this  arrangement 
is  calculated. 

To  facilitate  the  following  algorithms,  we  also  need  to  make  certain  additional  assump- 
tions concerning  the  function  <t>.  Roughly  speaking,  these  assumptions  require  that  the  result- 
ing data  structure  can  facilitate  especially  fast  calculation  of  <j)(x)  for  points  x  which  either  lie 
above  all  the  lines  containing  the  segmenU  in  G,  or  lie  below  all  these  lines.  See  Condition  A 
below  for  a  more  precise  statement. 
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2.2.   The  Preprocessing-and-Query  Version 

We  begin  by  describing  the  data  structure  to  be  used.  Our  description  follows  closely 
the  one  given  in  [EGSh],  which  in  turn  is  an  adaptation  of  the  partition-tree  approach  to 
half -plane  range  searching  [HW]  and  its  extensions  in  [DE].  The  reader  is  referred  to  these 
papers  for  more  information  concerning  the  construction. 

Specifically,  let  /,  be  the  line  containing  e,,  for  i  =  l,  .  .  .  ,n.  Dualize  these  lines  (in  the 
manner  used  in  [HW],  [EGSh]),  to  obtain  a  collection  G*  of  points  {/*,...  ,/*}.  Each  query 
point  p  then  becomes  a  dual  query  line  p  * .  We  construct  a  partition-tree  T  for  G  *  as  follows. 
With  each  node  v  of  T  we  associate  a  convex  subset  Cv  oi  the  dual  plane,  and  the  subset  G* 
of  the  dual  points  lying  in  Q^.  Let  n^  denote  the  cardinality  of  Gj.  If  nv=  1.  then  v  is  a  leaf 
of  T.  Otherwise  the  construction  of  T  below  v  continues  as  follows.  Take  the  set  G*  and 
apply  to  it  some  partitioning  scheme,  consisting  of  a  convex  subdivision  H^  of  Qy  having  cer- 
tain useful  properties  for  the  desired  point  location.  We  will  use  here  the  scheme  of  [HW] 
which  is  based  on  e-nets,  and  which  yields  the  second-best  query  performance  known  to  date. 
(The  best  is  a  recent  technique  of  Chazelle  and  Welzl  [We],  [CW].  This  technique  does  not 
seem  at  present  amenable  to  the  modifications  we  need  to  introduce  for  our  applications,  but 
has  been  successfully  applied  in  [Ag]  using  different  methods.)  Specifically,  for  some  arbi- 
trary fixed  6>0,  the  partition  Hy  =  Hi  '\%  &  convex  planar  subdivision  of  Q^  having 
A/  =  0((—  log  — )"*)  faces,  such  that  any  line  /  cuts  at  most  c  =  0((—  log  — )^)  faces  of  H^, 

€  €  €  € 

and  the  total  number  of  points  of  Gt  in  the  faces  cut  by  /  is  at  most  en.  Such  a  subdivision 
always  exists,  and  can  be  constructed  in  constant  randomized  time  [HW].  We  then  create  M 
children  wj,  .  .  .  ,wm  of  v,  one  per  face  of  H^,  associate  with  each  child  Wj  the  corresponding 
;'-th  face  Qy^j  of  H^,  and  the  subset  G'^  of  the  points  of  Gt  in  Q^j,  and  continue  this  partition- 
ing process  recursively  at  each  Wj,  thereby  obtaining  the  desired  tree  T. 

The  tree  T  is  used  to  process  a  query  point  p,  or  rather  its  dual  line  p  ',  as  follows.  We 

propagate  p*    through  T,   starting  at  the  root  of  T.   For  each  node  v  of  T,   let  ^"[p)  = 

2    <t>f(p)-    Suppose  p*  reaches  some  node  v  of  r  during  propagation.    If  v  is  a  leaf,  then 

^"{p)  is  simply  4),(/j),  where  e  is  the  single  element  of  G^  Otherwise,  by  the  properties  of 
Hy  mentioned  above,  p*  cuts  at  most  c  faces  of  H^,  and  "misses"  all  the  other  faces.  Let  Wj 
be  a  child  of  v.  Up*  cuts  the  face  Q^j  of  H^  associated  with  Wj  (this  will  be  reflected  by  say- 
ing that  p  *  "reaches"  Wj),  then  we  obtain  <t)*^(p)  recursively,  by  continuing  the  propagation 
of  p  *  past  Wj.  Up*  misses  Q^j,  we  need  an  ad-hoc  procedure,  to  be  discussed  below,  for 
obtaining  that  value.    Then  we  obtain  ^^(p)  by  straightforward  summation 

Vip)  =    2  <*>"^(P)  . 
>-» 

Note  that  when  p  '  misses  the  face  Qy^  of  a  node  w,  then  the  point  p  must  lie  either  above  all 
the  lines  /,  containing  the  segments  e,  in  G^.,  or  below  all  these  lines.  Thus  we  will  assume 
that,  using  appropriate  additional  preprocessing,  there  is  a  fast  way  to  compute  ^"^{p)  for 
such  points  p.  This  special-purpose  procedure  will  vary  from  one  application  to  another,  but 
in  all  applications  described  below  we  will  be  able  to  carry  out  this  computation  in  O(log  n^,) 
time  (as  noted  below,  we  can  also  tolerate  here  a  much  slower  computation).  Moreover,  the 
additional  preprocessing  needed  at  w  to  facilitate  this  fast  computation  will  usually  be  either 
0(n,^)  or  0(n,.,  log  n,,,),  and  will  consume  only  linear  space.  To  be  able  to  state  our  results  in 
general  but  precise  terms,  we  formalize  these  conditions  as  follows. 
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Condition  A:  Given  a  set  G  of  n  segments  in  the  plane,  one  can  preprocess  it  in  time 
0(n  log'^n)  (for  some  small  fixed  integer  r)  into  a  linear-size  data  structure,  so  that,  for 
any  query  point  p  lying  either  above  all  the  lines  containing  the  segments  in  G,  or  below 
all  these  lines,  one  can  calculate  <{>(/')  in  O(log  n)  time. 

Assuming  Condition  A,  we  can  easily  bound  the  query  time  Q  (n)  as  follows  (see  [HW] 
for  more  details). 

G(l)  =  0(1). 

G(n)  =  O (log /i)  +   2  CK)  .      '«>1. 
i  =  l 

c 

where  2  "'  ^  *"• 

i-l 

We  cl£iim  that  Q{n)  =  0{n'*)  for  any  y>2l2>.  This  is  proven  inductively,  as  in  [HW], 
using  the  inequality 

and  the  fact  that  log  n  =  0{n'^).  Since  •y>2/3  and  c  =  0((—  log  — )^),  it  follows  that  if  c  is 

chosen  sufficiently  small,  Q  {n)  =  O(n^)  will  satisfy  the  required  inequality.  (Note  that  the 
same  argument  still  applies  if  we  replace  log  n  by  any  function  which  is  0(n'').  This  allows  us 
to  relax  Condition  A  accordingly,  if  necessary.) 

As  to  the  preprocessing  time  and  space  of  the  algorithm,  it  is  easily  checked  that  the 
space  used  by  the  data  structure  is  0(n  log  n)  (the  partition  tree  T  itself  has  only  linear  size, 
and  the  auxiliary  structures  at  each  node  of  T,  as  given  by  Condition  A,  require  total  addi- 
tional 0(n  log  n)  storage).  Similarly,  the  time  needed  to  construct  the  tree  and  the  auxiliary 
structures  is  easily  seen  to  be  0(n  log'^'^'n).  Note  however  that  since  we  use  the  e-net  con- 
struction of  [HW],  processing  of  each  node  of  T  also  involves  a  (constant- time)  randomized 
selection  of  an  appropriate  c-net,  which,  with  high  probability,  produces  a  set  with  the 
desired  properties.  See  [EGSh]  for  more  details  concerning  this  issue,  and  also  for  an  alter- 
native fully  deterministic  construction  of  a  similar  partition  tree,  following  the  construction  in 
[EW],  which  yields  only  slightly  inferior  query  time. 

Moreover,  at  the  expense  of  slightly  complicating  the  query  processing,  we  can  reduce 
the  space  to  0(n),  with  a  similar  saving  of  0(log  n)  in  the  preprocessing  time,  following  the 
technique  of  [DE].  SPecifically,  we  store  the  auxiliary  data  structures  only  every  a  log  n  lev- 
els. To  compute  4)(p)  at  a  node  v  that  p  misses,  we  explore  a  complete  subtree  rooted  at  v 
and  having  depth  at  most  a  log  n  until  we  reach  a  level  where  the  auxiliary  data  structures 
are  stored.  We  then  compute  4)(p)  at  each  node  of  the  subtree  at  that  level  (obviously,  p 
misses  those  nodes  as  well),  and  sum  up  the  results  to  obtain  4)(p)  at  v.  The  cost  of  this 
extended  operation  is  O(log  n-2'''°«'')=0(n*),  for  any  8>0,  provided  a  is  chosen  appropri- 
ately small.  A  slight  modification  of  the  above  analysis  shows  that  the  query  time  is  still 
O(n^),  for  any  ■y>0,  and  the  space  requirement  is  now  only  0(n),  because  the  auxiliary 
structures  are  stored  at  only  a  constant  number  of  levels.  For  the  same  reason,  the  prepro- 
cessing time  also  decreases  by  a  factor  of  0(log  n). 

In  summary,  we  have  shown 
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Theorem  2.1:  Given  a  collection  of  n  segments  in  the  plane,  and  an  associated  function  <j>  as 
above  which  satisfies  Condition  A,  one  can  preprocess  these  segments  in  (randomized)  time 
0(n  log'n)  into  a  data  structure  of  size  0{n),  so  that,  given  any  query  point  p,  one  can  calcu- 
late <j>(p)  in  time  0{n^),  for  any  ■y>2/3. 

2.3.   The  Batched  Version 

We  now  turn  to  consider  the  batched  version  of  our  implicit  point  location  problem. 
The  set-up  is  the  same  as  assumed  in  the  preceding  subsection,  except  that  now  we  know  in 
advance  the  m  points  pi,  .  .  .  ,Pm  at  which  we  want  to  calculate  the  function  <j>.  We  use  the 
approach  of  [EGSh],  in  which  the  partition  tree  T  is  constructed  in  a  "customized"  style, 
catering  only  to  the  query  lines  p* ,  .  .  .  .pj,. 

Specifically,  we  construct  T  as  follows.  As  above,  let  /,  denote  the  line  containing  «,, 
for  each  i  =  l,  .  .  .  ,n.  We  construct  recursively  a  partition  tree  T  each  of  whose  nodes  v 
represents  a  convex  polygonal  region  Q^  of  the  dual  plane,  the  corresponding  subset  Gy  of 
the  rtv  segments  «,  for  which  the  dual  points  /*  lie  in  Q^  (we  denote  this  set  of  dual  points  by 
G* ),  and  the  set  Ly  of  the  m^  points  p;  whose  dual  lines  p*  cut  Q^  (we  denote  the  set  of  these 
dual  lines  by  L*)-  If  m^  <  nj,  we  expand  the  tree  T  at  v  as  follows.  We  choose  a  random 
sample  R  ol  r  oi  the  lines  in  L*  (where  r  is  some  constant  parameter  to  be  chosen  later),  and 
draw  and  triangulate  the  arrangement  Ay  of  these  sample  lines;  we  consider  only  the  portion 
of  Ay  lying  within  Cv  For  each  of  the  M  =  0{r^)  faces  /  of  this  portion  of  A^,  we  create  a 
child  w  of  V,  take  Q^  to  be  /,  and  obtain  G^,  Ly^  according  to  the  above  definition.  The  €-net 
theory  of  [HW]  is  easily  seen  to  imply  that,  with  high  probability,  the  number  m^  of  lines  p  * 

fHy  _ 

in  Lh.  is  at  most  0( log  r)  for  all  children  w.    If  my  a  ni  we  do  not  continue  the  construc- 

r 

tion  of  T  below  w,  pass  back  to  the  primal  plane,  construct  there  the  arrangement  A  {G^)  of 
the  Hy  segments  «,  in  Gy,  and  assign  to  each  face  (or  subface)  of  A  (Gy)  the  constant  value  of 
(J)"  over  that  face.  It  is  then  straightforward  to  collect  the  faces  of  A  (Gy)  that  contain  the  m^ 
points  Pi  that  have  reached  v,  and  obtain  the  value  of  4>^  at  each  of  these  points.  All  this  can 
be  easily  done  in  time  0{nl  +  Wy  log  n^)  =  0{m^  log  /ly),  and  space  0{m^). 

For  each  child  w  of  v,  we  collect  all  the  points  p,(.Ly  whose  dual  lines  have  missed  the 
face  Cw.  and  place  them  in  a  special  "bucket"  B^.  associated  with  w.  As  above,  if  p ,,  .  .  .  ,p, 
are  in  B^  then,  for  each  i,  the  point  p,  lies  either  above  all  the  lines  Ij  containing  the  seg- 
ments ej  €  Gy,  or  below  all  these  lines.  Condition  A  implies  that  the  values  <i>^ip),  forp  iBy^, 
can  all  be  calculated  in  total  time  0{by^  log  n^),  where  by^  =   |5h,|. 

For  an  inner  node  v  of  T,  we  obtain  (^"(pj)  for  all  the  points  Pj(.L„,  by  propagating 
these  lines  to  the  children  w,,  .  .  .  ,wi^  of  v,  obtaining  <|>"''(p;)  either  recursively,  if  p^€Lh,,, 
or  from  the  bucket  By^,  otherwise.  Then  the  values  of  4>^  at  all  these  points  is  readily  obtained 
by  summing  these  partial  results,  as  above. 

It  follows  that  the  time  T{m,n)  needed  to  calculate  the  value  of  <J),  for  a  given  collection 
of  n  segments,  at  m  given  points,  satisfies  the  following  recurrence  (in  which  we  assume  that 
the  random  sampling  at  each  node  of  T  does  indeed  produce  an  e-net  of  the  desired  type, 
which  will  be  the  case  with  high  probability): 

T{m,n)  =  0(m  log  n)  ,       if/nSn^ 

M 
T{m,n)  =    2  ^("'(.'•i)  +  0(m  log  n  -t-  n  log''n)  ,      otherwise 
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where  m,,  n,  satisfy  the  following  conditions. 
u 

(a)  2  "'  =  "• 

1  =  1 

(b)  m,  =  0( —  log  r)  for  each  i. 

Hence,  as  in  [EGSh],  we  obtain 

Theorem  2.2:  The  batched  version  of  the  implicit  point  location  problem,  as  formulated 
above,  can  be  solved  in  (randomized)  time 

0(m2'3-S„2'3  +  2S   +   ;„  log  „   +    „  log"--*-!/!) 

for  any  8>0. 

Remark:  As  already  noted  in  [EGSh],  Theorems  2.1  and  2.2  together  show  the  power  of 
prior  knowledge  of  the  points  p,  to  be  "located".  Indeed,  when  m  is  reasonably  large,  the 
above  procedure  spends  roughly  0((n^/m)"^^*  +  log  n)  time  per  point,  which  tends  to 
O(log  n)  as  m  approaches  n^. 

2.4.   An  example  -  point-triangle  containment  queries 

Let  A],  .  .  .  ,A„  be  n  (possibly  intersecting)  triangles  in  the  plane.  We  wish  to  prepro- 
cess  them  so  that,  given  any  query  point  x,  we  can  determine  quickly  whether  x  lies  in  the 

union  U   A,  of  these  triangles. 

This  containment  problem  fits  nicely  into  the  implicit  point  location  mold  presented 
above,  as  follows.  Let  A  be  a  triangle;  for  each  edge  e  of  A  let  B  («)  denote  the  semi-infinite 
trapezoidal  strip  of  points  lying  directly  below  e.  Define  a  function  4),  in  the  plane  so  that 
<|>,»«0  outside  B{e),  and  <j>,=  -t-l  on  B{e)  if  A  lies  below  the  Une  containing  e,  otherwise 
<!>,■- 1  on  B(e). 

It  is  easily  checked  that  if  «],e2>«3  ^^^  the  edges  of  A  then 


In  other  words,  if  we  put  <j>(x)  =  "^  <t>f(jt),  where  the  summation  extends  over  all  the  edges 


x€A 
otherwise 


of  the  triangles  A,,  then  4)(x)>0  if  and  only  if  x  is  contained  in   U   A,.  As  a  matter  of  fact, 

1  =  1 

<{)(x)  gives  the  number  of  triangles  A,  containing  x.  We  can  therefore  reformulate  our  prob- 
lem as  that  of  calculating  4)(p)  at  a  query  point  p,  and  testing  whether  <j>(p)>0.  We  thus 
want  to  apply  the  previous  analysis  to  the  collection  G  of  edges  «],...  ,«3„  of  the  triangles 

A,. 

To  do  so,  we  need  to  show,  for  any  given  subcoUection  Gv  of  G,  of  cardinality  n^,  that 

(i)      (j)"  (as  defined  in  Section  2.2)   is  readily  obtainable  during  calculation  of  the  entire 
arrangement  of  the  segments  in  G^; 

(ii)     4)^  satisfies  Condition  A. 

To  establish  the  first  property,  we  calculate  a  modified  arrangement  formed  by  the  n^ 
segments  in  G„  and  of  the  2n^  vertical  rays  extending  downwards  from  their  endpoints.  Cal- 
culation of  this  arrangement  still  takes  0(nj)  time,  and  by  definition  <j)  is  constant  on  each  of 
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its  faces,  and  can  be  obtained,  in  constant  time  per  face,  as  the  arrangement  is  being  con- 
structed. 

As  to  the  second  property,  let  p  be  a  point  lying  above  all  the  lines  containing  the  seg- 
ments in  Gv;  then  ^^(p)  =  0  by  definition.  On  the  other  hand,  if  p  lies  below  all  these  lines, 
we  can  obtain  <|>'^Cp)  as  follows.  We  project  the  2n^  endpoints  of  the  edges  in  Gv  onto  the  x- 
axis  and  sort  them  in  increasing  order.    Now,  given  such  a  point  p,  it  is  easily  checked  that 

<*»''(p)  =  S{«;  :«;€G„pe«;} 

where  e*  is  the  jt-projection  of  tj  and  where  e^  is  the  non-zero  value  of  <(>,..  Note  that  the 
sum  on  the  right-haud  side  is  constant  over  each  interval  between  two  successive  projected 
endpoints,  and  we  can  obtain  these  values  of  <i>^  by  maintaining  a  running  value  Z  of  this  sum 
as  we  scan  the  list  of  projected  points  from  left  to  right.  Initially,  2  =  0.  When  we  scan  the 
left  endpoint  of  some  e* ,  we  add  e^  to  2,  and  when  we  scan  the  right  endpoint  of  e*  we  sub- 
tract €j  from  2.  We  store  the  sorted  list  of  projected  endpoints,  together  with  the  associated 
values  of  ^y.  Thus,  given  a  query  point  p  lying  below  all  the  lines  containing  the  segments  of 
Gv,  we  can  obtain  4)^(p)  by  a  binary  search  through  this  list.  This  establishes  Condition  A 
(with  preprocessing  time  0(_n^  log  n^)). 

Remark:  In  general,  the  three  edges  of  a  triangle  do  not  "stay  together"  as  their  dual  points 
are  partitioned  among  the  nodes  of  the  tree  T.  As  a  result,  the  values  of  the  "intermediate" 
functions  4>^  at  inner  nodes  v  of  T  need  not  bear  any  relationship  to  the  original  containment 
problem. 

We  thus  obtain 

Theorem  2.3:  (i)  The  preprocessing  version  of  the  point-triangle  containment  problem  for  a 
collection  of  n  triangles  in  the  plane  can  be  solved  in  0(n  log  n)  randomized  expected 
preprocessing  time,  0(n)  space,  and  0{n^)  query  time,  for  any  -y>2/3. 

(ii)  The  batched  version  of  this  problem  for  n  triangles  and  m  points,  can  be  solved  in  ran- 
domized expected  time 

0(m2'3-»H2/3+25  +  „  logn  +  f.  log^n) 

for  any  8>0. 

Remark:  By  straightforward  modifications  of  the  above  procedure,  we  can  obtain  similar 
procedures  for  counting  how  many  triangles  A,  contain  a  query  point  x  (this  number  is  simply 
«!>(')).  or,  assuming  each  triangle  is  assigned  some  weight  w,,  for  determining  the  sum  of  the 
weights  of  the  triangles  that  contain  a  query  point  x  (this  is  achieved  by  an  appropriate 
modification  of  the  edge  functions  <J>,),  etc. 

3.    Ray  Shooting  in  Non-simple  Polygons 

We  now  consider  another  major  application  of  the  techniques  of  Section  2.  Let  P  be  a 
(connected)  polygonal  region  having  n  edges.  We  wish  to  preprocess  P  so  that,  given  any 
query  ray  p  emerging  from  any  point  inside  P,  we  can  determine  quickly  the  first  intersection 
of  p  with  dP  (or  determine  that  p  does  not  intersect  dP,  which  might  be  the  case  when  P  is 
unbounded).  If  P  is  simply  connected,  there  exist  efficient  techniques  for  solving  this  prob- 
lem [CG],  [GHLST]  which  require  0(T(n))  preprocessing  time,  C>(n)  space,  and  have 
0(log  n)  query  time  (here  j{n)=0{n  log  log  n)  is  the  time  needed  to  triangulate  an  n-sided 
simple  polygon  [TV]).  However,  if  /»  is  not  simply  connected,  the  problem  appears  to  be 
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more  difficult.  For  example,  Suri  and  O'Rourke  [SO]  have  studied  the  restricted  problem  of 
ray  shooting  into  a  non  simply  connected  polygon  P  (or  actually  of  finding  the  visible  portion 
of  P)  from  a  given  edge  of  P.  They  show  that  the  portion  of  P  visible  from  an  edge  is  a 
polygonal  region  that  can  be  described  as  the  union  of  O(n^)  triangles,  which,  in  the  worst 
case,  can  have  n(n^)  edges  on  its  boundary.  We  will  show,  nevertheless,  that  the  analysis  of 
[SO]  is  over-pessimistic,  and  that  this  more  general  version  of  the  problem  can  be  attacked 
using  our  partition-tree  approach,  and  can  thus  be  solved  relatively  efficiently. 

For  technical  reasons,  and  without  loss  of  generality,  we  will  consider  only  rays  that  are 
rightward-directed.  Let  ej,  .  .  .  ,e„  be  the  edges  of  P.  We  use  a  multi-level  data  structure, 
such  that  each  level  is  based  on  a  partition-tree  for  a  different  type  of  data  points  and  query 
lines.  The  primary  level  is  a  partition-tree  T  for  the  set  G '  =  {/*,...,/*}  of  dual  points  to 
the  lines  /,  containing  the  edges  e,  of  P.  The  dual  q*  oi  the  origin  q  of  the  query  ray  p  will 
serve  as  a  query  line  to  be  propagated  through  T.  At  each  node  v  of  T  there  is  associated  a 
subset  Gv  of  the  edges  of  P  (so  that  the  set  G*  of  the  duals  of  the  lines  containing  the  seg- 
ments in  Gv,  is  the  set  of  points  /*  lying  in  the  convex  face  Gv  of  the  partitioning  that  v 
represents),  and  the  recursive  goal  at  v  is  to  find  the  nearest  intersection  of  p  (if  any)  with 
any  of  the  segments  in  Gv  If  v  is  a  leaf  of  T  this  task  is  trivially  done  in  constant  time, 
because  Gv  contains  just  one  segment.  If  v  is  an  inner  node  of  T,  then  the  answer  to  the 
query  at  v  is  easily  obtained  in  constant  time,  provided  we  already  know  the  answer  at  each 
child  of  V. 

There  are  two  ways  to  get  the  desired  solution  at  each  child  w  of  v.  If  the  query  line  q* 
reaches  w  (in  the  sense  that  it  intersects  the  corresponding  face  Cv)  then  p  and  q*  are  pro- 
pagated to  w  and  the  answer  is  obtained  by  recursive  processing  of  w.  Otherwise,  q* 
"misses"  w,  and  a  special  procedure  is  needed  to  obtsiin  the  answer  to  the  query  at  w. 

Note  that  if  9*  misses  w,  then  the  point  q  must  either  lie  above  all  the  lines  containing 
the  segments  in  Gh-,  or  lie  below  all  these  lines.  Assume  without  loss  of  generality  the  first 
case  (the  second  case  is  handled  in  a  completely  symmetric  fashion). 

To  handle  this  special  situation,  we  define  a  linear  ordering  <  on  the  segments  in  G^ 
which  is  consistent  with  our  ray-shooting  problem,  in  the  sense  that  if  our  ray  p  (emerging 
from  a  point  q  as  above)  hits  two  segments  *,*'  in  G^  so  that  the  intersection  of  p  with  e  lies 
to  the  left  of  its  intersection  with  e' ,  then  e<e' .  This  will  enable  us  to  apply  binary  search 
(in  a  non-trivial  way,  though)  to  find  the  first  intersection  of  p  with  these  segments. 

Similar  orderings  for  other  classes  of  planar  objects  have  been  investigated  by  Guibas 
and  Yao  [GY]  and  by  Rosenstiehl  [Ro].  Some  of  the  proof  techniques  used  below  are 
adapted  from  [GY]. 

3.1.    Linear  Ordering  of  Segments  for  Ray  Shooting 

Let  «!,...,«„  be  n  segments  in  the  plane  with  pairwise  disjoint  (relative)  interiors. 
For  each  e,,  its  top  side  e^  is  that  side  of  e,  visible  from  a  point  at  y  =  -t-oo,  and  its  bottom 
side  ef  is  the  opposite  side  (if  e,  is  vertical,  we  define  its  top  side  to  be  its  left  side,  and  take 
its  bottom  side  to  be  undefined). 

DeHnition:  A  (non-vertical)  line  /  is  said  to  hit  «,  at  its  top  side,  or  simply  to  hit  e,* ,  if  /  inter- 
sects e,  transversally  at  some  point  x,  and,  slightly  to  the  left  of  x,  I  lies  above  e,.  /  is  said  to 
hit  eT  when  a  symmetric  condition  holds.  (If  e,  is  vertical,  any  line  /  intersecting  e,  is  said  to 
hit«r.) 
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Definition:  e,  <j  tj  if  there  exists  a  (non-vertical)  line  /  hitting  both  ef  and  ef  such  that  its 
intersection  x,  with  «,  lies  to  the  left  of  its  interset:  on  Xj  with  ej,  and  such  that  /  does  not  hit 
any  bottom  side  of  any  other  tj^  at  a  point  between  x,  and  Xj. 

Definition:  «,  <v  tj  if  there  exists  a  vertical  line  intersecting  both  «,  and  tj  such  that  its  inter- 
section with  «,  lies  above  its  intersection  with  ej. 

Lemma  3.1:    <v   is  acyclic. 

Proof:  This  is  also  proved  in  [EGS],  by  showing  that  any  cycle  can  be  shortened  by  omitting 
the  segment  with  the  leftmost  left  endpoint.  We  will  show,  using  a  different  proof,  that  in 
any  finite  collection  E  of  segments,  there  always  exists  a  segment  that  is  minimal  under  <v  . 
This  segment  e  is  taken  to  be  the  one  whose  right  endpoint  (or  top  endpoint  if  it  is  vertical) 
can  "see  the  sky"  (i.e.  the  upward  vertical  ray  from  that  endpoint  does  not  hit  any  segment 
of  E),  and  which  is  such  that  its  right  endpoint  is  leftmost  among  all  segments  in  E  satisfying 
this  condition.  To  see  that  e  is  minimal,  suppose  to  the  contrary  that  there  exists  some 
e'  <v  «  in  E.  Then  it  is  easily  checked  that  the  right  endpoint  of  e'  lies  above  e.  Let  e"  be 
the  segment  whose  right  endpoint  lies  above  e  and  is  rightmost  among  all  segments  in  E  with 
this  property.  Then  it  is  easily  checked  that  the  right  endpoint  of  e"  can  see  the  sky  (and  lies 
to  the  left  of  the  right  endpoint  of  e),  a  contradiction  which  completes  the  proof.  D 

Definition:  «,  <,  «y  if  «,  and  ej  have  non-overlapping  x-projections  and  the  projection  of  e, 
lies  to  the  left  of  that  of  ej. 

Definition:  e,  <v*/x  «;  (this  is  our  main  order  relationship  which  will  also  be  denoted  shortly 
as  e,  <  ej)  if  either  e,  precedes  ej  in  the  transitive  closure  <v*  of  <v  ,  or  e,  and  ej  are  not 
related  by    <v»    and  «,  <^  ej. 

Lemma  3.2:    <v»/x   is  a  linear  order. 

Proof:  First  observe  that  every  pair  «,,  ej  are  related  by  this  relationship.  Indeed,  if  e,  and  ej 
are  not  related  by  <^*  ,  then  their  x-projections  are  non-overlapping,  which  implies  that 
they  are  related  by    <;,  . 

It  remains  to  be  shown  that  <y'u  is  acyclic.  Let  us  first  define  yet  another  relation- 
ship: «,  <(-v*)/vx  ej  if  e,  and  ej  are  not  related  by  <v*  and  e,  <^  ej\  thus  e,  <v*/i  «;  if  and 
only  if  either  e,  <„•  «;  or  «/  <(-v*)/\r  <;•  By  definition,  <v»  is  transitive:  e  <^.»  e'  and 
e'  <v*  e"  implies  e  <y»  e" .  We  claim  that  <(-v*)/u  is  also  transitive.  Suppose 
«  <(.v*)Ax  «'  and  e'  <(.v»)ax  «"•  Clearly  e  <^  e'  and  e'  <^  e" ,  and  thus  e  <^  e" .  The  claim 
will  thus  follow  if  we  show  that  e  and  e"  are  not  related  under  <v*  .  Suppose  to  the  con- 
trary that  e  <y»  e" .  Thus  there  is  a  chain  of  relationships  «  =  e,o  <v  «,i  <v  •  "  •  <v  <n 
<v  en  +  i  =  «".  Without  loss  of  generality  we  can  assume  this  chain  to  be  maximal,  i.e.  that 
there  is  no  longer  chain  of  such  inequalities  between  e  and  e" .  Since  e,  <v  e,  ^j  in  this  maxi- 
mal chain,  there  exists  a  vertical  segment  6^  connecting  a  point  Op  on  e,p  to  a  point  bp  on 
'ip  +  v  ^^^^  "p  lying  above  bp,  so  that  Sp  does  not  intersect  any  segment  of  our  collection. 
Consider  the  path  y  obtained  by  the  concatenation 

y  =  8oll*oail|8,||fcia2||  •  •  •  ||fc,_,aj|8,  . 

-y  is  a  simple  path  that  starts  at  a  point  to  the  left  of  e' ,  ends  at  a  point  to  the  right  of  e' ,  and 
does  not  intersect  <•'.  Therefore  y  must  pass  either  above  e'  or  below  e' .  But  by  definition  of 
<v  it  is  easily  checked  that  in  the  first  case  we  have  e  <v»  e' ,  and  in  the  second  case 
«'  <v'  «".  a  contradiction.  A  similar  contradiction  is  obtained  by  a  completely  symmetric 
argument  in  the  case  e"  <v«  e. 
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Now  we  can  prove  the  acyclicity  of   <v»/x  ■    Suppose  to  the  contrary  that  there  is  a  cycle 

<ii  <  «,2  <  •  •  •  <  «.„  <  e,i 

and  that  this  cycle  is  of  minimal  length  m.  First  of  all,  it  is  easily  seen  that  m  =3.  Indeed, 
m  =  2  is  ruled  out  by  definition,  and  any  cycle  of  length  >3  can  be  shortcut.  Thus  we  can 
take  our  cycle  to  be  «  <  e'  <  e"  <  e.  Now  each  of  the  three  relationships  along  the  cycle  is 
either  <v*  or  <(.v*)ai  •  Not  all  three  relationships  can  be  <v*  ,  because  that  would  imply 
that  <v  is  cyclic,  contrary  to  Lemma  3.1.  Also,  not  all  three  can  be  <(,v*)ai  .  because  that 
would  imply  in  particular  e  <x  e'  <x  e"  <,  e,  which  is  clearly  impossible.  Suppose  two  of 
the  links  are  <v*  and  one  is  <(,v*)ax  •  Without  loss  of  generality  we  can  write 
e  <v*  e'  <v*  e"  <(.v*)ai  «•  But  by  transitivity  we  obtain  e  <v«  e"  which  contradicts,  by 
definition,  the  relationship  e"  <(_v*)/u  «•  Similarly,  if  e  <(_v*)ai  «'  <(-v*)Ar  «"  <v*  «.  then 
by  transitivity  of   <(_v*)ftx   proved  above,  we  again  reach  a  contradiction 

This  shows  that   <v*/x   is  total  and  acyclic,  and  thus  a  linear  order.  □ 

Lemma  3.3:  If  «  <,  «'  then  e  <v*/j  «'• 

Proof:  By  definition,  there  exists  a  line  /  hitting  both  e  and  e'  on  their  top  side,  so  that  its 
intersection  point  x  with  e  is  to  the  left  of  its  intersection  x'  with  e' ,  and  such  that  /  hits  no 
bottom  side  between  x  and  x' .  Without  loss  of  generality,  we  can  also  assume  that  /  hits  no 
top  side  of  a  segment  between  x  and  x' .  Indeed,  if  /  hits  the  top  sides  of  e^j,  .  .  .  ,eji^ 
between  x  and  x'  in  this  order  from  left  to  right,  then  it  suffices  to  prove  our  lemma  for  each 
of  the  pairs  e<ejy,ej^<ej2,  .  .  .  ,eji^<e',  and  then  obtain  our  claim  from  the  transitivity  of 

<v*/x   • 

If  e  and  e'  have  overlapping  x-projections,  then  they  are  related  by  <v  and  it  is  easy  to 
check  that  e  <v  «',  which  implies  e  <v»/i  e' .  Otherwise  we  must  have  e  <x  e' .  If  e  and  e' 
are  not  related  by  <v*  ,  or  if  e  <v*  «',  then  again  e  <v*/x  «'•  Thus  we  can  assume  e'  <v»  e, 
i.e.  we  have  a  chain  of  the  form 

where  again  we  assume  this  chain  of  inequalities  to  be  maximal. 

Let  tr  denote  the  unbounded  polygonal  curve  obtained  by  the  concatenation  of  the 
downward-directed  vertical  ray  r  emerging  from  the  left  endpoint  a  of  e,  the  segments  ax, 
xx' ,  x'b' ,  where  b'  is  the  right  endpoint  of  e' ,  and  the  upward-directed  vertical  ray  r'  emerg- 
ing from  b' .  IT  partitions  the  plane  into  two  regions  R^,  L„,  lying  respectively  to  its  right  and 
to  its  left  (as  we  traverse  ir  from  r  to  r'). 

We  claim  that  none  of  the  segments  e,j,  .  .  .  ,e,^  intersects  it.  Indeed,  none  of  them 
can  intersect  ax,  xx' ,  or  x'b' .  If  some  segment  e,p  intersects  r,  then  by  definition  we  have 
e  <v  tip,  implying  a  cycle  in  <v  and  contradicting  Lemma  3.1.  A  symmetric  argument 
implies  that  no  segment  can  intersect  r' . 

As  in  the  proof  of  Lemma  3.2,  let  7  be  the  simple  directed  polygonal  path  obtained  by 
following  the  (maximal)  chain  of  inequalities  e'  <y  tj^  <v  •  •  •  <v  <•  Since  the  first  link  of 
■y  lies  below  e'  (thus  in  /?„)  and  the  last  link  of  -y  lies  above  e  (thus  in  L„),  it  follows  that  -y 
must  intersect  -n .  As  argued  above,  none  of  the  segments  «,_  can  intersect  it.  Consider  the 
first  intersection  z,  encountered  along  -y,  of  this  curve  with  it.  Then  z  lies  along  one  of  the 
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vertical  subsegments  of  -y,  connecting  a  point  a€e,  with  a  point  b^tj^^  for  some  p,  so  that  a 
lies  above  b.  But  the  only  portion  of  it  that  ab  can  intersect  is  the  segment  xx' ,  and,  by 
assumption,  a  lies  in  Rj,,  hence  below  xx' .  Thus  b  also  lies  below  xx' ,  so  the  intersection 
between  ab  and  xx'  is  impossible,  a  contradiction  which  completes  the  proof  of  the  lemma,  o 

Corollary:    <,   is  a  partial  order,  and    <v*/i   is  a  linear  extension  of  it. 

We  now  turn  to  describe  an  algorithm  for  calculating  <v»/x  .  The  idea  is  quite  simple. 
Using  a  vertical  line-sweep  algorithm,  it  is  easy  to  obtain  all  0{n)  pairs  of  segments  («,«') 
related  by  <v  ,  such  that  there  exists  a  line  /  meeting  both  e  and  e'  and  the  portion  of 
between  these  segments  docs  not  meet  any  other  segment.  Clearly,  the  collection  E  of  such 
pairs  is  a  subset  of  <v  whose  transitive  closure  coincides  with  <v*  •  We  next  apply  a  topo- 
logical sorting  procedure  to  £,  so  as  to  obtain  a  linearization  of  <v*  ,  and  adapt  it  so  that 
the  linearization  it  produces  is  actually  <v*/x  •  The  standard  topological  sort  procedure  [Kn] 
maintains  a  subset  /  of  the  given  segments,  all  of  whose  predecessors  (if  any)  have  already 
been  placed  in  the  output  linear  order,  and  then  selects  at  each  step  an  element  of  /  as  the 
next  element  in  the  linear  order.  Note  that  each  pair  of  segments  in  /  have  non-overlapping 
x-projections.  We  modify  the  procedure  so  that  it  maintains  the  set  /  as  a  priority  queue 
(ordered  by  x-coordinate),  and  always  selects  from  it  the  smallest  element  as  the  next  ele- 
ment to  be  placed  on  the  output  list. 

It  is  easy  to  show,  using  induction  on  the  order  in  which  the  segments  are  being  placed 
on  the  output  list,  that  the  resulting  order  coincides  with  <v*ix  ■  To  see  the  general  induc- 
tion step,  consider  the  current  value  of  /,  and  let  e  be  the  smallest  segment  (in  the  x-ordering) 
in  /.  It  suffices  to  show  that  e  is  the  smallest  element,  in  the  order  <v»/j  ,  among  all  the  still 
unprocessed  segments.  Suppose  the  contrary.  Then  there  exists  another  unprocessed  seg- 
ment *'  such  that  «'  <v*/x  <■  Without  loss  of  generality  we  can  assume  that  *'€/,  because 
every  unprocessed  segment  must  have  a  predecessor,  in  the  order  <v*  ,  in  /.  But  then 
e  <x  e' ,  so  that  we  must  have  «'  <v»  e.  But  since  e'  has  not  yet  been  processed,  e  cannot 
belong  to  /,  a  contradiction  which  proves  our  claim. 

The  time  complexity  of  the  algorithm  is  clearly  0{n  log  n),  which  arises  from  the  cost 
of  the  initial  sweeping  algorithm,  and  from  the  maintenance  of  the  priority  queue  /  in  the 
topological  sorting.  All  other  steps  take  linear  time.  We  thus  have  in  summary 

Theorem  3.4:  Given  a  collection  of  n  segments  in  the  plane  with  disjoint  interiors,  one  can 
compute,  in  time  0{n  log  n),  a  hnear  order  of  them  which  extends  their  "rightward  shooting 
order"    <,    as  defined  above. 

3.2.   The  Ray  Shooting  Procedure,  Continued 

Let  us  now  return  to  the  special  case  of  our  shooting  problem.  The  ordering  <  has  been 
defined  above  in  more  generality  than  we  actually  need  for  this  special  case,  because  the  rays 
that  we  consider  now  cannot  hit  the  bottom  side  of  any  segment  in  Gw  This  fact  yields  the 
following  simple  corollary  which  is  essential  to  our  procedure. 

Lemma  3.5:  Suppose  <,e'€G»,  and  e  <,  *'.  If  a  rightward-directed  ray  p  emerges  from  a 
point  lying  above  all  the  lines  containing  the  segments  in  G^  and  hits  e,  then  either  p  does  not 
hit  e'  or  it  hits  it  at  a  point  lying  to  the  right  of  the  intersection  of  p  with  e. 
Proof:  If  p  were  to  intersect  e'  to  the  left  of  its  intersection  with  e,  then  since  both  these 
intersections  must  hit  the  top  sides  of  e,e' ,  and  since  p  cannot  hit  the  bottom  side  of  any  seg- 
ment  in   Gh,,    it   follows   by    definition    that   we   would    have   e'  <,  e,    a   contradiction    that 
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completes  the  proof,  n 

We  will  use  this  result  to  obtain  a  divide-and-conquer  strategy,  based  on  repeated  appli- 
cations of  the  following  blocking  detection  procedure: 

We  are  given  n  segments  «],  .  .  .  ,e„  and  we  want  to  preprocess  them  so  that,  given  any 
query  (rightward-directed)  ray  p  originating  at  a  point  lying  above  all  the  lines  contain- 
ing the  given  segments,  we  wish  to  determine  quickly  whether  p  intersects  any  of  the 
segments  «i,  .  .  .  ,e„. 

The  solution  to  this  problem  is  a  straightforward  application  of  the  implicit  point  loca- 
tion technique  described  in  Section  2.  Specifically,  let  /  denote  the  line  containing  the  query 
ray  p.  It  is  easily  verified  that  p  intersects  a  given  segment  «,  if  and  only  if  the  line  /  hits  the 
top  side  of  «,.  Moreover,  passing  to  the  dual  plane,  this  condition  is  equivalent  to  the  condi- 
tion that  the  dual  point  /*  lies  in  the  right  wedge  W,  formed  by  the  two  dual  lines  to  the  end- 
points  of  «,.  Thus  our  blocking  detection  query  is  equivalent  to  testing  whether  /*  is  con- 
tained in  at  least  one  of  these  wedges,  and,  as  we  have  seen  in  Theorem  2.3,  this  can  be 
determined  in  time  O^n"^)  (for  any  'y>2/3),  after  0(n  log  n)  randomized  expected  prepro- 
cessing and  using  O  (n  log  n)  space. 

Having  this  solution  at  hand,  we  solve  our  shooting  problem  as  follows.  Let 
ei  <  ei  <  •  '  ■  <  e„  (where  n  =  n„)  be  the  linear  ordering  <v*/x  of  ttie  segments  in  Gh-. 
The  preprocessing  for  this  problem  consists  of  constructing  the  above  data  structure  for  the 
blocking  detection  problem  for  the  first  half  {«i,  .  .  .  ,c„/2}  of  the  segments  in  Gh.,  followed 

by  recursive  preprocessing  of  each  of  the  two  halves  {«i 'n/i)'  {«n/2  +  i>  •  ■  ■  ,'n}  of  the 

set  Gh..    Clearly,  this  preprocessing  requires  0(n  log^n)  time  and  uses  0(n  log  n)  space. 

However,  we  can  apply  once  again  a  space-saving  trick  akin  to  that  of  [DE].  That  is,  we 
split  («],  .  .  .  ,e„)  into  n*  consecutive  subsequences,  each  of  size  n'~*,  apply  the  preprocess- 
ing step  of  the  blocking  detection  procedure  for  each  subsequence  separately  (other  than  the 
last  one),  and  then  preprocess  recursively  each  subsequence  for  the  shooting  problem.  More- 
over, in  each  recursive  application  we  split  the  corresponding  subsequence  into  n*  subse- 
quences, where  n  is  the  size  of  the  original  set.  In  this  way,  there  is  only  a  constant  number 
of  recursive  levels.  It  is  easily  verified  that  the  overall  cost  of  this  recursive  preprocessing  is 
still  0(n  log  n),  and  the  total  space  is  0{n). 

To  process  a  shooting  query  involving  a  (rightward-directed)  ray  p  emerging  from  a 
point  lying  above  all  the  lines  containing  the  segments  of  Gy^,  we  first  test,  in  0{n^)  time, 
whether  p  is  blocked  by  the  first  subsequence  of  G^.  If  so,  then  Lemma  3.5  implies  that  we 
can  ignore  in  this  case  the  remaining  subsequences  of  G^,  and  we  thus  continue  recursively  to 
query  the  first  subsequence.  Otherwise,  we  apply  the  blocking  detection  procedure  to  the 
second  subsequence  of  Gh,,  and  continue  this  way  until  we  find  a  subsequence  that  blocks  p, 
and  then  continue  the  query  recursively  within  that  subsequence.  At  the  end  of  this  process, 
we  reach  a  subset  of  Gh  consisting  of  a  single  segment,  which  we  return  as  the  output  to  our 
shooting  query.    The  overall  query  time  is 

0(n^**  +  (Y)^**  +  (j)^*'+  •  •  •  )  =  0(n>**) 

which  can  also  be  written  as  0{n^)  for  another,  yet  still  arbitrarily  small,  ■y>2/3. 

Our  final  stroke  is  to  apply  the  space-saving  trick  of  [DE]  once  again,  by  preprocessing 
the  sets  Gh,  just  once  every  h  levels,  where  h  is  chosen  appropriately  as  above.  The  above 
arguments  imply 
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Theorem  3.6:  Given  a  polygonal  region  P  with  n  edges,  one  can  preprocess  it  in  time 
0{n  log  n)  into  a  data  structure  of  size  0(n),  using  which  one  can  answer  ray  shooting 
queries  in  time  OCn""),  for  any  y>2/3. 

3.3.    Ray  Shooting  in  an  Arbitrary  Arrangement  of  Segments 

The  preceding  techniques  can  be  generalized  to  handle  the  more  difficult  case  of  ray 
shooting  in  aa  arrangement  of  an  arbitrary  collection  of  (possibly  intersecting)  line  segments. 
Let  {«],...,«„}  be  such  a  collection.  We  prepare  a  data  structure  similar  to  the  one  used 
above.  At  the  top  level  we  use  the  same  partition  tree  T  as  above.  In  answering  a  shooting 
query,  the  processing  of  leaves  of  T,  and  the  merging  of  query  outputs  at  children  of  a  node  v 
of  r  to  obtain  the  output  at  v  itself,  are  the  same  as  above.  The  part  that  needs  modification 
is  in  the  processing  of  a  node  w  of  T  that  the  dual  line  q*  (where  q  is  the  origin  of  the  query 
ray)  has  missed. 

Assume  without  loss  of  generality  that  q  lies  above  all  the  lines  containing  the  segments 
in  Cw-  Let  R  denote  the  intersection  of  the  upper  half-planes  bounded  below  by  these  lines. 
Clearly  q^R  and  the  interior  of  R  is  disjoint  from  any  segment  in  Gm>.  Let  F  denote  the  face 
of  the  arrangement  of  the  segments  in  G,^  which  contains  R  (if  all  segments  in  G^  are 
bounded,  F  is  actually  the  (unique)  unbounded  face  of  the  arrangement).  As  shown  in  [PSS], 
[EGSh],  the  boundary  of  F  consists  of  0(n,^a(n^^))  subsegments,  and  can  be  calculated  in 
time  O  (n^a(ny,.)  log^n^).  Once  F  is  available,  the  remainder  of  the  construction  proceeds 
exactly  as  above,  except  that  this  time  it  is  applied  to  the  0(n^a(n^))  subsegments  (which 
have  pairwise  disjoint  relative  interiors)  forming  the  boundary  of  F. 

Embedding  the  cost  of  this  additional  preprocessing  into  that  of  the  entire  algorithm,  we 
easily  obtain 

Theorem  3.7:  Given  a  collection  of  n  segments  in  the  plane,  one  can  preprocess  it  in  time 
0(na(n)  log  n)  into  a  data  structure  of  size  0(n  a(n)),  using  which  one  can  answer  ray 
shooting  queries  (in  which  we  seek,  for  a  given  query  ray,  the  first  segment  it  hits)  in  time 
O(n^),  for  any  -y>2/3. 

4.    Implicit  Hidden  Sarface  Removal  Algorithms 

As  our  next  application,  we  study  certain  cases  of  the  hidden  surface  removal  problem, 
which  can  be  stated  as  follows.  Given  a  collection  of  objects  in  3-dimensional  space,  and  a 
viewing  point  a,  we  wish  to  calculate  the  scene  obtained  by  viewing  these  objects  from  a.  In 
other  words,  we  want  to  produce  some  representation  of  this  scene,  in  which  we  can  deter- 
mine, for  each  ray  emerging  from  a,  the  first  surface  of  an  object  hit  by  that  ray  (or  deter- 
mine that  the  ray  does  not  his  any  object).  For  the  sake  of  simplicity,  we  will  assume  that  the 
given  objects  are  all  polyhedral,  whose  boundaries  consist  of  n  triangular  faces  altogether, 
and  that  they  do  not  intersect  one  another.  This  problem  can  be  easily  formulated  in  the  con- 
text of  implicit  point  location,  as  follows.  The  given  triangular  faces  are  projected  towards  a 
onto  a  certain  "plane  of  view",  and  the  arrangement  of  their  projections  is  calculated.  To  each 
face  of  this  arrangement  one  then  assigns  the  triangle  nearest  to  a  whose  projection  contains 
that  face,  and  this  augmented  arrangement  constitutes  the  desired  solution  to  the  hidden  sur- 
face removal  problem.  To  obtain  an  actual  image  for  display  purposes,  one  can  then  take 
each  pixel  in  the  plane  of  view,  and  locate  the  face  of  the  arrangement  containing  that  pixel, 
thereby  obtaining  the  object  seen  at  that  pixel.  Known  solutions  to  this  problem  run  in  time 
and  space  0{n^)  [De],  [MK],  which  is  of  course  worst-case  optimal  if  exphcit  calculation  of 
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the  entire  image  is  desired.  We  will  use  the  implicit  algorithms  of  Section  2  to  obtain 
improved  solutions.  This  is  not  the  first  instance  of  an  implicit  solution  to  this  problem.  It 
was  shown  recently  in  [CS]  that  if  the  objects  to  be  viewed  form  a  polyhedral  terrain  (i.e.  a 
surface  cut  by  each  verticeil  line  in  exactly  one  point)  then  we  can  obtain  an  implicit  represen- 
tation of  its  image  in  preprocessing  time  and  space  0(na(n)logn)  (where  a(n)  is  the 
extremely  slowly  growing  inverse  Ackermann's  function),  and  query  time  O(log  n). 

We  will  consider  the  following  special  case  of  the  hidden  surface  removal  problem.  Let 
T  =  {Ai,  .  .  .  ,A„}  be  a  collection  of  n  horizontal  triangles  in  3-space  such  that  A,  lies  in  the 
plane  z=i  (or,  more  generally.  A,  lies  in  a  plane  z  =  a,  where  0<ai<  •  •  •  <a„).  The  prob- 
lem is  to  determine,  for  each  point  p  in  the  xy-plane,  the  first  (i.e.  lowest)  triangle  Ay  hit  by 
the  vertical  line  passing  through  p  (or  to  determine  that  no  such  triangle  exists).  What  distin- 
guishes this  special  case  from  the  general  surface  removal  problem  is  that  in  this  case  there 
exists  a  natural  total  ordering  (by  their  z-coordinates)  of  the  triangles  in  T,  with  the  property 
that  if  a  vertical  line  hits  both  A,  and  A^  with  A,  lying  below  Ay,  then  A/  <  Ay.  Our  results 
will  actually  apply  to  any  instance  of  the  hidden  surface  removal  problem  in  which  such  an 
ordering  exists. 

We  also  consider  the  following  simpler  problem.  Given  the  same  set-up  as  above,  we 
want  to  determine,  for  any  point  p  in  the  ry-plane,  whether  p  is  blocked  by  the  triangles 
Ai A„,  meaning  that  the  vertical  line  through  p  intersects  at  least  one  of  these  trian- 
gles. In  other  words,  we  W£mt  to  determine  whether  p  lies  in  the  union  of  the  projections 
A*,  .  .  .  ,AJ  of  the  given  triangles  onto  the  j:y-pleme.  For  this  blocking  problem,  no  order 
among  the  given  triangles  need  be  assumed.  In  fact,  it  applies  to  an  arbitrary  collection  of 
(even  intersecting)  triangles  in  3-space.  By  Theorem  2.3,  we  immediately  obtain 

Theorem  4.1:  (i)  The  preprocessing  version  of  the  visibility  blocking  problem  for  a  collection 
of  n  triangles  in  3-space  can  be  solved  in  0{n  log  n)  preprocessing  time,  0{n)  space,  and 
Oin"^)  query  time,  for  any  ■y>2/3. 

(ii)  The  batched  version  of  this  problem  for  n  triangles  and  m  points,  can  be  solved  in  time 

0(m2'3-*n2'3-^28  4-  „  log  n  -(-  n  log^n) 

for  any  8>0. 

Remark:  We  leave  it  as  an  open  problem  whether  the  results  of  Section  2  can  be  applied 
directly  to  solve  the  original  hidden  surface  elimination  problem.  The  difficulty  we  face  is 
that  while  membership  of  a  query  point  p  in  a  projected  triangle  is  obtained  by  adding  edge 
functions,  finding  the  lowest  triangle  above  p  involves  calculating  the  minimum  of  the  indices 
of  the  triangles  lying  above  p.  This  mixture  of  addition  and  minimum  operations  does  not 
seem  to  be  amenable  to  the  techniques  of  Section  2,  which  tend  to  "scramble"  the  order  in 
which  these  operations  are  applied. 

To  solve  our  original  hidden  surface  removal  problem,  we  continue  in  much  the  same 
way  as  in  Section  3.  Consider  first  the  preprocessing  version.  We  split  the  set  T  of  triangles 
into  two  subsets  T,,  Tj,  so  that  T,  =  {Ai,  .  .  .  .A^/j}  contains  the  lower  half  of  the  triangles, 
and  T2  =  {A„,2+i.  •  •  •  ,A„}  contains  the  upper  half  of  the  triangles.  We  apply  the  blocking 
detection  preprocessing  algorithm  to  the  triangles  in  T,,  and  then  proceed  recursively  with 
preprocessing  T]  and  T2  separately.  To  process  a  query  point  p,  we  first  test  whether  T, 
blocks  p.  If  so,  we  continue  to  process  the  query  within  the  structure  of  Tj;  otherwise  we 
continue  with  the  structure  of  T2.  Thus  in  log  n  such  steps  we  find  the  lowest  trjangle  block- 
ing p   (if  one  exists).  This  approach  uses  0{n  log^n)   preprocessing  time  and   0(n  log  n) 
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space,  and  has  query  time  0(n"*)  for  any  7>2/3.  Again,  the  trick  of  [DE]  can  be  used  to 
reduce  the  preprocessing  time  and  space  to  0(n  log  n)  and  0{n),  respectively. 

The  batched  version  is  handled  similarly.  We  split  the  set  T  into  two  subsets  Tj,  T2  as 
above,  and  then  apply  the  batched  version  of  the  blocking  detection  procedure  to  the  m  given 
points  and  to  Tj.  Suppose  m  1  of  the  points  are  found  to  be  blocked  by  T],  and  nij  =  m  -  nix 
are  not  blocked.  We  then  calculate  recursively  the  lowest  triangle  in  Tj  above  each  of  the 
first  m  1  points,  and  similsirly  calculate  the  lowest  triangle  in  T2  above  each  of  the  remaining 
m2  points.  If  during  the  recursion  the  number  of  points  or  the  number  of  triangles  becomes 
sufficiently  small  (e.g.  ^1),  then  we  solve  the  problem  by  brute  force,  in  maximum  time 
0(m+n). 

Hence,  if  T{m,n)  denotes  the  maximum  time  to  solve  our  problem  for  m  points  and  n 
triangles,  we  have 


T{m,n)  =  0{m  +n),       if  m^l  or  nsl 


r(m,n)  =      max 

m  \+  m2^rn 


7-(m,,|)  +  r(m2,|) 


+  C(m,n),     otherwise 


where 


Q(m,n)  =  0((m2'3-8„2'3  +  28  +  m  +  n  log  n)  log  n) 


r(m,n)   =    0((m2'3-»„2/3  +  28    +   ^  ^^^2^    ^    „  ,^g3„) 


for  any  6>0. 
Proposition  4.2: 

for  any  8>0. 

Proof:  It  is  clear  that  the  contribution  of  the  second  and  third  terms  of  Q(m,n)  to  T(m,n)  is 
bounded  by  the  second  and  third  terms  of  this  expression.  As  to  the  first  term,  the  inequality 
is  an  immediate  consequence  of 

„2/3-5(i)2/3.26   +    ^2/3-»  (|)2/3.28   ^   („  ^    +   ^^)2/3-8„2.3.26  .  ^^ 

<   ^2/3-5^2/3  +  28 


Theorem  4.3:  (i)  The  preprocessing  version  of  the  hidden  surface  removal  problem  for  an 
ordered  collection  of  n  triangles  in  3-space  can  be  solved  in  0(n  log  n)  randomized  expected 
preprocessing  time,  0(n)  space,  and  O(n^)  query  time,  for  any  7>2/3. 

(ii)  The  batched  version  of  this  problem  for  n  triangles  and  m  points,  can  be  solved  in  ran- 
domized expected  time 

0(m2'3-»„2/3  +  25   +   „   ,^,g2„    +    „  jog3„) 

for  any  8>0. 

Remark:  Recently  and  independently,  MuUer  [Mu],  [SML]  has  studied  other  variants  of  the 
ray  shooting  problem  in  three  dimensions.  In  particular,  he  has  looked  at  cases  where  the 
given  objects  are  all  iso-oriented  horizontal  rectangles,  but  the  query  ray  can  be  in  arbitrary 
directions.  He  has  applied  somewhat  simUar  partitioning  techniques,  and  obtained  similar, 
though  somewhat  worse,  performance  bounds.  This  is  due  to  the  following  two  reasons,  both 
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of  which  can  be  rectified  using  our  approach,  (i)  He  has  used  the  conjugation  tree  technique 
of  [EW]  rather  than  the  €-net  approach  of  [HW];  and  (ii)  In  the  batched  version,  he  has  used 
standard,  non-customized  trees. 

5.    Polygon  Placement  Queries 

Consider  next  the  following  final  application  of  our  techniques.  Let  P  be  any  polygonal 
object  (not  necessarily  simply  connected)  having  k  sides,  and  let  A],  ...  ,A„  be  a  collection 
of  (possibly  intersecting)  triangles  in  the  plane.  We  wish  to  preprocess  the  data  so  that,  given 
any  query  placement  of  (a  fixed  reference  point  within)  P ,  we  can  quickly  determine  whether 
at  this  placement  P  hits  any  of  the  triangles,  or  is  otherwise  "free"  (we  assume  that  P  is  kept 
at  its  initial  orientation).  Problems  of  this  sort  arise  in  a  variety  of  contexts  (see  Chazelle 
[Ch]).  If  P  is  convex  (and  the  "obstacles"  do  not  intersect  each  other),  then  there  exist  effi- 
cient solutions  to  the  problem  ([CD],  [Fo],  [LS],  [KLPS],  [BZ]),  but  the  general  case  has  no 
comparably  efficient  solutions.  The  general  approach  to  this  problem  is  to  form  the  Min- 
kowski differences 

Ki  =  ^i-P  =  {x-y.   x€A„  yiP) 

(where  x—y  denotes  vector  subtraction),  and  then  calculate  their  union 

n 

K  =    \J   Ki  . 

i=l 

Assuming  the  reference  point  on  the  given  initial  placement  of  P  is  the  origin,  then  a  place- 
ment of  P  with  that  reference  point  at  some  point  z  (we  say  in  short,  placement  of  P  at  z)  is 
free,  if  and  only  \i  ziK. 

Now  if  we  triangulate  P  into  0(t)  triangles,  say  /'),...  ,Pt.,  we  can  represent  K  as  the 
union  of  0{kn)  Minkowski  differences  of  the  form  Kij=Ai  —  Pj,  each  of  which  is  a  convex  /- 
gon,  for  some  /s6.  Thus  our  blocking  detection  procedure  can  be  applied  to  these  Kjj,  to 
yield  the  following  result. 

Theorem  5.1:  With  preprocessing  time  0(kn  log  kn)  and  space  0(kn),  one  can  obtain  a  data- 
structure,  using  which  one  can  determine,  for  a  given  query  placement  of  P,  whether  at  this 
placement  P  is  free  or  not,  in  time  0{{kny)  for  any  "y>2/3. 

Remark:  We  can  also  obtain  efficient  solutions  to  the  batched  version  of  the  problem,  as 
above. 

How  significant  is  this  result?  If  the  obstacle  triangles  A,  intersect  one  another,  it  seems 
that  our  implicit  approach  is  superior  to  any  explicit  alternative  technique.  However,  sup- 
pose the  triangles  A,  have  pairwise  disjoint  interiors.  Then,  if  P  is  convex,  the  previously 
mentioned  results  [CD],  [Fo],  [LS],  [KLPS],  [BZ]  show  that  the  complexity  of  K  is  only 
0{kn),  and  that  it  can  be  calculated  in  a  variety  of  efficient  techniques,  the  best  being  of  time 
complexity  0(kn  log  kn)  [LS].  Using  standard  point  location  techniques,  we  can  then  process 
each  placement  query  in  time  0(log  kn). 

Now,  this  observation  yields  an  alternative  approach  to  the  case  of  a  general  P. 
Namely,  we  decompose  P  into  triangles  P,,  .  .  .  ,/'t,  and  then  apply  the  above  techniques  to 
each  triangle  P,  separately.  Given  a  query  placement  z  of  P,  we  locate  z  in  each  of  the  * 
maps  obtained  for  the  individual  parts  P,  of  P.  Then  z  is  free  if  and  only  if  it  is  free  in  each 
of  these  maps.  This  gives  us  preprocessing  time  0(kn  log  n),  space  0(kn),  and  query  time 
O(tlogn). 
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Comparing  this  query  time  to  our  solution,  it  is  easily  checked  that  our  method  becomes 
superior  when  k>n^. 

Remark:  Both  approaches,  that  of  the  k  separate  queries  and  our  implicit  point  location 
approach,  also  apply  to  cases  where  P  is  defined  as  the  union  of  k  arbitrary  triangles  or  line 
segments.  Thus  even  in  cases  of  this  kind  where  the  "explicit"  complexity  of  P  is  O(i^)  (e.g. 
P  is  a.  kxk  grid  of  horizontal  and  vertical  line  segments),  we  can  determine  whether  a  place- 
ment of  P  is  free  or  not  within  the  same  bounds  noted  above. 

6.    Conclusion 

The  main  purpose  of  this  paper  has  been  to  point  out  the  versatility  and  rather  wide 
applicability  of  the  partition  tree  technique.  We  have  demonstrated  this  by  developing  a  gen- 
eral technique  for  "implicit"  point  location  in  arrangements  of  overlapping  polygonal  objects, 
and  by  three  other  applications  to  ray  shooting  and  containment  problems.  We  believe  that 
there  are  many  additional  potential  applications  of  this  method. 
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